Hilbert-Schmidt Geometry of n-Level Jakobczyk-Siennicki 
Two-Dimensional Quantum Systems 

Paul B. SlateJ] 

ISBER, University of California, 
Santa Barbara, CA 93106 

(Dated: February 1, 2008) 

Abstract 

Jakobczyk and Siennicki studied iiyodimensional sections of a set of (generalized) Bloch vec- 

o : 

■ tors corresponding to n x n density matrices of two-qubit systems (that is, the case n = 4). 



in 
o 
o 

o 

Q 



o 



They found essentially five different types of (nontrivial) separability regimes. We compute the 
Euclidean/Hilbert-Schmidt (HS) separability probabilities assigned to these regimes, and conduct 



> 

(N 

' parallel iwo-dimensional sectional analyses for the higher-level cases n = 6, 8, 9 and 10. Making 

o ■ 

. use of the newly-introduced capability for integration over implicitly defined regions of version 



5.1 of Mathematica — as we have also fruitfully done in the n = 4 i/iree-parameter entropy- 



maximization-based study quant-ph/0507203 — we obtain a wide-ranging variety of exact HS- 



■ probabilities. For n > 6, the probabilities are those of having a partial positive transpose (PPT). 

^ : 

Q^, For the n = 6 case, we also obtain biseparability probabilities; in the n = 8, 9 instances, bi-PPT 

probabilities; and for n = 8, tri-PPT probabilities. By far, the most frequently recorded probabil- 
' ity for n > 4 is ^ ~ 0.785398. We also conduct a number of related analyses, pertaining to the 

(one-dimensional) boundaries (both exterior and interior) of the separability and PPT domains, 
and attempt (with quite limited success) some exact calculations pertaining to the 9-dimensional 
(real) and i5-dimensional (complex) convex sets of two-qubit density matrices — for which exact 
HS-separability probabilities have been conjectured, but not yet verified. 
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I. INTRODUCTION 

There has been considerable recent interest 0, 0, Q| in understanding how one can, from 
the spherical coordinate point-of view, generalize to n-level quantum systems {n > 2) the 
familiar Bloch ball representation of the two-level quantum systems (n = 2) — in which the 
pure states form the bounding spherical surface ("Bloch sphere") of the unit ball in three- 
dimensional Euclidean space. Kimura and Kossakowski have expressed the generalized Bloch 
representation of an n x n density matrix in the form eq. (3)] 

P=^/n + ^s£r'(trpA.)A„ (1) 

where /„ is the identity operator, and the Aj's are the (n^ — 1) orthogonal generators of 
SU{n), forming a basis of the set of all the linear operators with respect to the Hilbert- 
Schmidt inner product. 

An interesting application of these concepts was made by Jakobczyk and Siennicki (JS) 
^. They examined all those two-qubit (n = 4) systems describable as iwo-dimensional 
sections of sets of (generalized) Bloch (coherence ^) vectors. (The totality of 4 x 4 density 
matrices, on the other hand, comprises a /z/iteen-dimensional convex set — the nxn density 
matrices being [n^ — l)-dimensional in nature — so thirteen of the fifteen SU{4) orthogonal 
generators [Gell-mann matrices] are assigned null weight in the JS r?, = 4 analyses. That is, 
thirteen of the fifteen coefficients, (trAj) in the expansion are zero.) 

Since there were only two parameters involved in each of their scenarios, JS were able 
to present planar diagrams depicting the feasible regions, as well as those subsets of these 
regions composed of separable states. In their Fig. 1, JS exhibited thirteen possible types 
of parameter domains. Further, in their Fig. 2, they showed six different (nontrivial) sepa- 
rability scenarios (two of which — labelled "EF)" and "FE)" by JS - are simply geometric 
reflections of one another). 

We will, firstly (sec. |n|), in this study, evaluate the sizes (areas) of these six (two- 
dimensional) domains and nontrivial subdomains, in terms of the Hilbert- Schmidt (HS) 
metric — a task JS did not exphcitly address. (The HS-distance between two density 
operators pi,p2 is defined as ^y^T{pl — p'l)^ P, eq. (2.3)].) Then, we extend the JS analyses 
to the cases r?, = 6 (sec. IIIIIl . 8 (sec. IIV|1 . 9 (sec. EI) and 10 (sec. IVIIl . in which various 
mM/i2partite — as opposed to simply bipartite scenarios can arise. For all these instances. 
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except n = 10, we additionally obtain the HS-lengths of the (one-dimensional) boundary 
states (that is, those with degenerate spectra) and the HS-probabilities that states lying 
on this boundary are separable. (Motivated by our extensive numerical results given in 
0] and ^, Szarek, Bengtsson and Zyczkowski have recently proved "that the probability 
to find a random state to be separable equals 2 times the probability to find a random 
boundary state to be separable, provided the random states are generated uniformly with 
respect to the Hilbert-Schmidt (Euclidean) distance" [s^). Also, we compute in certain cases, 
the HS-lengths of the (interior) boundaries dividing one domain of interest from another. 
(The interior states generically have nondegenerate spectra.) Then (sec. IVIIIl . we undertake 
some analyses involving three (rather than two) parameters. These prove to be much more 
problematical in nature (cf. P]). 

We also report, at the end fsec. lVIin]) . some initial steps in an attempt to determine exact 
upper bounds for the HS- volumes of the separable 5-dimensional real and i5-dimensional 
complex 4x4 density matrices. (Only in these computations — in order to compare our 
formulas with known HS-volumes of separable and nonseparable steps — do we not take 
the HS- volume element to be unity.) 

Our computations in this paper were greatly facilitated by a new feature of the pro- 
gramming language Mathematica (version 5.1) — the capacity to integrate over implicitly 
defined regions. (This feature was also employed by us in j^, in a somewhat related two- 
qubit context, in which the Jaynes maximum-entropy principle was employed.) We, first, 
found explicit forms for the n eigenvalues of the various n x n matrices and for their par- 
tial transposes. Then we required, in the several integrations, using the new feature (thus, 
saving us from the laborious task of having to specify large numbers of particular integra- 
tion limits and do corresponding detailed bookkeeping), simply that these eigenvalues be 
nonnegative. This ensured that we either had, in fact, the requisite density matrices and/or 
positive partial transposes (PPT) of density matrices. 



II. THE QUBIT-QUBIT CASE n = 4 OF JAK0BCZYK AND SIENNICKI 

For the (geometrically-refiected) scenarios that JS labeled "EF)" and "FE)", we have 
found (Table H]) that the Hilbert-Schmidt volume (cf. j^) of separable and nonseparable 
states is and of the separable states alone is |. So the corresponding separability 
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probability (taking ratios) is — elegantly — "75 ~ 0.707107. 

For the scenario "CK)", possessing a triangular separability domain, the total volume is 
and the separability probability is j^{9 + 2V^rr) ^ 0.828450. For "GH')", the total 
volume is ^-y/fTr and the separability probability, — (2^/2) ^ 0.825312. For "KG)", 
the total volume is and the separability probability is (the smallest) | + 2^ ~ 0.608998. 
For "HG')", the total HS-volume is and the HS-separability probability is (the largest 
of the five) ''"^'^I^X''^'^ 0.842035. 

TABLE I: 

JS scenario HS total vol. HS separable vol. HS sep. prob. num. approx. 

EF) and EE) ^ 1 73 0.707107 

CK) ^ ^ (3^/6 + 2V2tt) 0.828450 

GH') Ayj, ^(52V3 + 27y6sin-i(2^)) 0.825312 

HG') ^ ^(52^3 + 27^/6 sec-i (3)) "'""'"l^X"'^'^ 0.842035 

KG) ^ (3^/6 + 2^/27r) k + ^ 0.608998 

Now, let us present again most of these results (Table H} in the form of the array 0. We 
do so because we will also present all the results of our subsequent analyses below (for n > 4) 
in this manner (which we have found to be the most convenient for directly incorporating 
our large-scale Mathematica computer-generated analyses into this report). 

In the first column of Q are given the identifying numbers of a pair of Gell-Mann matrices 
(generators of SU{4)) — which, in fact, can be seen to fully agree with the numbering (and 
associated scenario-labelling) of JS y, d. 389] . (Here and further, we will always adhere to 
the conventional/standard numbering 1^ sec. HI] of the Lie generators of SU{n), so that 
our results should be reproducible/verifiable to others. We list the pairs in lexographic order, 
using the first pair as the representative for its equivalence class.) In the second column of (j21) 
are shown the number of distinct unordered pairs of SU (4) generators which share the same 
total (separable and nonseparable) HS volume, as well as the same separable HS volume, 
and consequently, identical HS separability probabilities. The third column gives us these 
HS total volumes, the fourth column, the HS separability probabilities and the last (fifth) 
column, numerical approximations to the exact probabilities (which, of course, we see — 
being probabilities — do not exceed the value 1). (Due to space/page width constraints, we 



were unable to generally present in these data arrays the HS separable volumes too, though 

they can, of course, be deduced from the total volume and the separability probability.) 

[ {3,6} A ^ 0.707107 ^ 

{6,8} 2 A^l^ 26V2-,2nan-n2V2) ^ ^25312 

{6, 15} 2 ^ ^ (9 + 2V37t) 0.828450 

{8,9} 2 ^ ''''''^X'^'^ 0.842035 

^{9,15} 2 ^ l + il 0.608998 



(2) 



Thus, twelve of the 210 = 15-14 possible unordered pairs of Gell-Mann matrices are asso- 
ciated with nontrivial (< 1) separability probabilities [4, p. 389]. 

A. Boundary states 



For the scenario associated with the pair of Gell-Mann matrices {3,6}, the HS-length 
of the (one-dimensional) boundary states (that is, those with degenerate spectra) is |, and 
of the bounding states which are separable, |. For the pair {6,8}, the analogous results 
are ^ ^ 1.06066 and 1; for {6, 15}, they are | and |; and for {8,9}, ^ and for a 
separability probability of boundary states of |. For the last {9, 15} of the five scenarios, 
we have and 1. Let us now present these results in the following array form (which we 
will adopt for our more extensive results further below): 



1 {3,6} 


3 
2 


1 
2 


1 
3 


0.333333 \ 


{6,8} 


3 


1 


3 


0.942809 


{6,15} 


2 
3 


1 
2 


3 
4 


0.75 


{8,9} 


3^/3 
4 


V3 
2 


2 

3 


0.666667 


V{9,15} 


2 

V3 


1 


V3 
2 


0.866025 / 



(3) 



B. Length of Separability-Nonseparability Interior Boundary 



In the following array, we present the HS-length of the common border separating the 
nonseparable (entangled) states from the separable ones. The states lying along this interior 
border generically have nondegenerate spectra. 

{3,6} {6,8} {6,15} {8,9} {9,15} 
i 1 1 



4 



(4) 
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III. THE QUBIT-QUTRIT CASE n = 6 



A. 3x2 Decomposition 



Moving on from the n = 4 case specifically studied by Jacobczyk and Siennicki to n = 6 
(cf. [g^), we compute the partial transposes of the 6x6 density matrices, corresponding to 
two-dimensional sections of the set of Bloch vectors. We, first, transpose in place the (2^) 
four 3x3 blocks of the density matrices. By the Peres- Horodecki criterion, such density 
matrices with positive partial transposes must be separable. 

We obtained the following results, presented in the same manner as 0. 



{1,13} 
{3,11} 
{3,13} 
{3,25} 
{8,13} 
{8,25} 



48 I 
8V2 



4 
4 
4 
4 
4 



{11,15} 4 

{11,24} 2 

{13,24} 2 

{13,35} 4 

{15,16} 4 

{16,24} 2 

{20,24} 2 

{24,25} 2 

{24,27} 2 

^{25,35} 4 



27 

4 
9 

27 

2 
3 

V2 
3 

4v^7T 

27 
25 



4 

1 

V2 
5 
6 

5 

4v/2 
1 

3 

1 I 3vl 

3 ~'~ 47r 



72 



25 



72 

Vi 

5 

32v^ 
81 



25 
144 



25 
144 

25 
27V2 

25 
27 

5 



(-2 + 5^5) 
^(5 + 3v^csc-i 
;^(9v^ + 4vr) 



4+5 sin 



TT 



TT 



, 5 / 



5tt 




92+75 sin-l 


1) 


7571 




^ 5^/5 
92+75 cos-i 


(!) 


80^5 
V5+3csc-i( 




37r 



0.785398 ^ 
0.707107 
0.833333 
0.883883 
0.577350 
0.816497 
0.746830 
0.863648 
0.979236 
0.886838 
0.879838 
0.549815 
0.685627 
0.821115 
0.903076 

0.504975 y 



(5) 



1. Boundary states 



In the following array, we list, first the scenario pair, then, the HS-length of the boundary 
states (those with degenerate spectra), then, the HS-length of those boundary states which 
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are separable, then, the separabihty probabihty and a numerical approximation to it. 



{1,13} 


2 

3 


U 


U 


u. 


To -1 -1 T 

{3,11} 


1 


1 
3 


1 
3 


0.333333 


To 1 Ol 

{3, 13} 


2 
3 


1 
3 


1 
2 


U.O 


{3,25} 


1 


1 

3 


1 

3 


0.333333 


{8, 13} 


2 

73 


1 

2^3 


1 
4 


0.25 


{8, 25} 


%/3 
2 


1 
V3 


2 
3 


O.oooooT 


111 1 c^l 


4 

3v/3 


2 

3 


V3 
2 


U.oDDUZO 




5 

12 


1 
3 


4 
5 


U.O 


{13,24} 


5 
6 


V5 
3 


2 

7E 


0.894427 


{13,35} 


2 
5 


1 

3 


5 
6 


0.833333 


il5 161 


3 


V 3 


3 
4 


0.75 


ilfi 941 


5 
6 


2 
3 


4 
5 


U.O 


{20,24} 


5 
6 

6 


2 

3 


4 
5 


0.8 


{24,25} 


vl 

6 


1 

5 


0.2 


{24,27} 


6 


3 


4 
5 


0.8 


1^{25,35} 


2 

V5 


2 
3 


3 


0.745356 y 
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2. Length of Separability-Nonseparability Interior Boundary 

In the following arrays, we present the HS-length of the common border separating the 
nonseparable (entangled) states from the separable ones for each specific scenario. 



{1,13} {3,11} {3,13} {3,25} {8,13} {8,25} {11,15} {11,24} 



1 1 

VS 2a/3 



(7) 



{13,24} {13,35} {15,16} {16,24} {20,24} {24,25} {24,27} {25,35} 

-\/f 1 2 Vl 1 



V5 2 
3 3 



B. 2x3 Decomposition 



Here, we compute the partial transposes of the same collection of 6 x 6 density matrices, 
corresponding to two-dimensional sections of the set of Bloch vectors, by transposing in 
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place the (3^) nine 2x2 blocks of the density matrices 
blocks as previously fsec. IIIl A]) . 
We obtained the following results. 



rather than the four (2^)3 x 3 



{3,6} 8 ^ ^ 



{4, 18} 
{6,8} 
{6,15} 
{8,9} 
{8,22} 
{8,29} 
{9,15} 
{15,22} 2 
{15,29} 2 
{18,24} 4 
{18,35} 4 
{24,25} 4 
\^ {25,35} 4 



48 I 



IT 



16 



27 

V2 
3 

1 

V3 

2 

3 

4V27r 



27 



32 



27 

32^2 
81 

25 
144 



In 



26^2+27 tan-l(2V2) 
27^? 

^ (9 + 2VSn) 

52+27^2 sec-i (3) 
48^6 

8 
9 

4 
3v^ 

1 _L V3 
3 ~'~ 27r 

1 
2 

VA 
2 

92+75 sin-i(l) 
75^? 



^ ^(5 + 3v/5csc-(^ 

92+75 cos- il'^'l 



25 
27\/2 

5 



80^5 
V5+3csc-i(^) 



0.707107 
0.785398 
0.825312 

0.828450 
0.842035 
0.888889 
0.769800 
0.608998 
0.500000 
0.866025 
0.685627 
0.886838 
0.903076 
0.504975 



There are now only fourteen rows, while in the preceding qubit-qutrit analysis (0) there were 
sixteen. In both analyses, though, there are 48 unordered pairs of Lie generators which yield 
HS separability probabilities equal to |. 
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1. Boundary states 

Here, we again present the results, restricting consideration to the boundary (degenerate 
spectra) states, in the same form as previously 0- 



1 {3,6} 


1 


1 
3 


1 

3 


0.333333 \ 


{4, 18} 


2 
3 








0. 


{6,8} 


1 

V2 


2 

3 


2v^ 
3 


0.942809 


{6,15} 


4 
9 


1 
3 


3 
4 


0.75 


{8,9} 


%/3 
2 


1 

V3 


2 

3 


0.666667 


{8,22} 


x/3 
2 


1 


2 
3 


0.666667 


{8,29} 


2 

V3 


1 

2^3 


1 
4 


0.25 


{9,15} 


4 
3v^ 


2 

3 


V3 
2 


0.866025 


{15,22} 




3 


1 

5 


0.2 


{15,29} 


3 


/? 

V ^ 


3 
4 


0.75 


{18,24} 


5 
6 


2 
3 


4 
5 


0.8 


{18,35} 


2 

5 


1 
3 


5 
6 


0.833333 


{24,25} 


6 


3 


4 
5 


0.8 


U25,35} 


2 

75 


2 
3 


V5 
3 


0.745356 / 
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2. Length of Separability-Nonseparability Interior Boundary 

In the following arrays, we present the HS-length of the common border separating the 
nonseparable (entangled) states from the separable ones for each specific scenario. 

I {3,6} {4,18} {6,8} {6,15} {8,9} {8,22} {8,29} 

I 1 2 2 2 _\_ _J_ J_ 

\ 3 3 3 3 2^3 2^3 

{9,15} {15,22} {15,29} {18,24} {18,35} {24,25} {25,35} 

1 fi y2 2 2 v5 1 

3 "^3 3 3 3 6 3 

C. Biseparable HS probabilities 

Now, we determine which of the two-dimensional set of 6 x 6 density matrices have positive 
partial transposes for both forms of partial transposition used in sec. IIII AIIII^ The results 
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obtained were: 





oo 


4 
9 


{3,6} 


12 


27 


{3,13} 


4 


4 

Q 

y 


{3,27} 


2 


Z / 


{6,8} 


2 


8 y 2 


16 151 


2 


27 




2 


%/2 
3 


-f8 13]- 


4 


2 
3 


|0, ZZj- 


9 


1 

V3 


l8 251- 


4 


V2 
3 


{8,29} 


2 


2 

Q 
O 


{9,15} 


2 


97 


{11,15} 


4 


97 


{11,24} 


2 


72 


{13,24} 


2 


72 


{13,35} 


8 


5 


{15,16} 


4 


32^2 
81 


{15,22} 


2 


27 


{15,29} 


2 


32V2 
81 


{16,24} 


2 


25 fb 
144 Y 2 


{18,24} 


6 


25 [5 
144 Y 2 


{24,25} 


2 


25 
27t/2 


{24,27} 


4 


25 
27\/2 


l^{25,35} 


6 


5 



TT 
4 


0.785398 




0.707107 


5 
6 


u. 000000 


5 

4^2 

26^2+27 tan^l (2^2) 
277r 


0.883883 
0.825312 


2^ (9 + 2V37r) 


0.828450 


52+27\/2sec~^(3) 
48^6 


0.842035 




0.577350 


9 


0.888889 


Vi 


0.816497 


4 

3a/3 


0.769800 


1 1 V3 

3 27r 


u.ouoyyo 


1 _|_ 3V3 


746830 


5 + 2Sin (gj 


U.oDoD4o 


(-2 + 5^/5) 


0.979236 


^(5 + 3v/5csc-^ (75). 


) 0.886838 


i (9\/3 + 47r) 


0.879838 


1 
2 


0.500000 


V3 
2 


0.866025 


4+5sm-l(|) 
57r 


549815 


92+75 sin-i(l) 

757r 


0.685627 


3(4+5cos-i(|)) 

16v^ 
92+75 cos-i(l) 

80a/5 
V5+3csc-i(^) 

37r 


0.724191 
0.903076 

0.504975 
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1. Boundary States 



Concerning the corresponding one-dimensional (exterior) boundary (degenerate spectra) 
states we found: 

^{1,13} I 0. ^ 

{3,6} 1 I I 0.333333 

{3,13} f I I 0.5 

{3,27} 1 i I 0.333333 

{6,8} -ij I ^ 0.942809 

{6,15} I i f 0.75 

{8,9} ^ 75 I 0.666667 

{8' 13} ^ ^ i 0.25 

{8,22} ^ ^ I 0.666667 

{8,25} ^ 75 I 0.666667 

{8' 29} ^ ^ i 0.25 

{9,15} i ^ 0.866025 

{11,15} ^ I ^ 0.866025 " ^^^^ 

{11,24} A 1 I 8 

{13,24} I ^ ^ 0.894427 

{13,35} III 0.833333 

(15, 16} 1^ ^ I 0.75 

{15,22} ^ ^ I 0.2 

{15,29} ^ I 0.75 

{16,24} I i t 0.8 

{18,24} I 1 t 0.8 

{24,25} ^ 0. 

{24,27} ^ 4^ t 0.8 

V {25, 35} ^ I f 0.745356 ) 
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2. Length of Biseparability-Nonbiseparability Interior Boundary 



In the following array, we present the HS-length of the common border separating the 
(generically nondegenerate) biseparable states from the non-biseparable ones for each specific 



scenario. 



{1,13} 


2 
3 


{3,6} 


2 
3 


{3, 13} 


2 
3 


{3, 27} 


2 

3 


{6,8} 


4 
3 


{b, 15| 


1 
i 


{8,9} 


V3 
2 


|8, 13} 


^/3 
2 




^/3 
2 




V3 
2 


18,29} 


V3 
2 


{9,15} 


1 


{11,15} 


4 
3 


{11,24} 


1 


{13,24} 


2v^ 

Q 
O 


{13,35} 


1 


{15,16} 




{15,22} 


3 


{15,29} 


3 


{16,24} 


1 


{18,24} 


4 

3 


{24,25} 


\/f 1 
6 


{24,27} 


6 


l^{25,35} 


2 
3 



10 
3 

\/lO 
3 



(13) 
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IV. THE CASE n = 8 



A. 4x2 Decomposition 

Here, we compute the partial transposes of the 8x8 density matrices, corresponding to 
two-dimensional sections of the set of Bloch vectors, by, first, transposing in place the (2^) 
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four 4x4 blocks of the density matrices. The results are 



( {1,20} 
{3, 18} 
{3,20} 
{3,36} 
{8,20} 
{8,22} 
{8,36} 
{8,42} 
{8,49} 
{8,53} 
{15,22} 
{15,49} 
{18,24} 
{18,35} 
{20,35} 
{20,48} 

{22,48} 
{22,63} 

{24,25} 
{25,35} 
{29,35} 
{35,36} 
{35,38} 
{36,48} 
{42,48} 
{48,49} 
{48,53} 
^{49,63} 



192 
4 



1 

4 

1 

3a/2 

1 

4 

1 

3 

8 

3V3 
16 

3 

3^3 
16 

3 
8a/2 

3 
8 

Vi 

3 

9 

25 
256 

9 



TT 



20 



20 



49 



192 

fVi 

192 

8 

25 
48\/2 

9V37r 
80 

9-v/37r 
80 

3 
5V2 

3 
5\/2 

49 
384 

49 
384 

49 
72 V2 

49 
72 a/2 

27r 
7 



1 

5 
6 

5 
4V2 
1 

^/3 
7 
9 

'J 

3 

9 

7 

3^6 

4 
3V3 

1 

2 

V3 
2 

92+75sin-i(|) 
75^? 

^ (|5 + 3v^csc-^ 

ik (-25 + 24^30) 
^ (^12 + 7V6sin-i 

Ti7 (-9 + 7V42) 
^ (|7 + 4v^csc-i 

92+75cos-i(|) 
80t/5 



3 

V5 



/ 2V1 



4 

V7 



\/5+3c 



37r 



14^5+27 CSC" 

i_ Vi 



24 



TT 



TT 



14\/5+27cos-i(|) 

We 

2y6+7sm-i(2^) 



77r 



22V6+49 



497r 



3(22^6+49 cos-l(l)) 

^+4csc-l(-^) 
47r 



0.785398 \ 
0.707107 
0.833333 
0.883883 
0.577350 
0.777778 
0.816497 
0.888889 
0.952579 
0.769800 
0.500000 
0.866025 
0.685627 
0.886838 
0.985680 
0.903278 

0.989529 
0.915544 

0.903076 
0.504975 
0.636783 
0.809819 
0.918793 
0.469522 
0.596820 
0.801610 
0.930129 
0.440596 J 
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1. Boundary states 



The results pertaining to the one-dimensional (exterior) boundary (generically degener- 
ate) states were: 



/ 


1 

2 








0. \ 


{3,18} 


3 
4 


1 

4 


1 

3 


0.333333 


{3,20} 


1 
2 


1 
4 


1 
2 


0.5 


{3,36} 


3 

4 


1 

4 


1 

3 


0.333333 


{8,20} 


V3 
2 


V3 
8 


1 
4 


0.25 


{8,22} 


ov "J 

8 


Vo 
8 


1 
i 

3 


0.333333 


{8,36} 


ov o 

8 


V <J 

4 


2 
3 


0.666667 


{8,42} 


3\/3 


V 

4 


2 
3 


0.666667 


{8,49} 


3V3 


\/3 
4 


2 

3 


0.666667 


{8,53} 


V3 
2 


8 


1 
4 


0.25 




5 
2V6 


1 

2^/6 


1 
5 


n 


/it; /IQ\ 


[2 
V 3 


2 


3 
4 


u. / 


{18,24} 


5 
8 


1 
2 


4 
5 


U.o 




3 

10 


1 
4 


5 
6 


U. 000000 




3 
5 








-f9n 4Rl 


7 

24 


1 

4 


6 
7 


u.ou / ±^0 


|22,48| 


7 
12 


2 




U.yzooz 




2 
7 


1 

4 


7 

8 


n eye: 
0.0/0 


{24,25} 


8 


Vi 
2 


4 
5 


U.o 


{25,35} 


3 

2^5 


1 

2 


V t> 
3 


0.745356 




3 
2^5 


1 
2 


V5 
3 


U. / 4000D 




3^ 
2 


4 


1 

6 


U.iDDDD/ 


{35,38} 


2 


4 


5 
6 


0.833333 


{36,48} 


7 
4\/6 


1 
2 


2V6 
7 


0.699854 


{42,48} 


7 

4^6 


1 

2 


2V6 
7 


0.699854 


{48,49} 


8 
8 


VI 

8 


1 
7 


0.142857 


{48, 53} 


4 


6 
7 


0.857143 


U49,63} 


2 


1 
2 


V7 
4 


0.661438 ) 



(15) 
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Here, we see the appearance of (fully entangled) domains that have no separable component, 
at all. 



B. 2x4 Decomposition 



Now, we compute the partial transposes of the 8x8 density matrices, corresponding to 
two-dimensional sections of the set of Bloch vectors, by transposing in place the (4^) sixteen 
2x2 blocks of the density matrices. We obtained the following results. 



{3,6} 12 
{4,18} 192 \ 
{6,8} 2 



1 

3V2 



{6, 15} 

{8,9} 
{8,22} 
{8,29} 
{9,15} 
{15,22} 4 
{15,29} 4 
{18,24} 4 
{18,35} 4 
{24,25} 4 
{24,46} 2 
{24,57} 2 

{25,35} 4 
{35,46} 2 
{35,57} 2 

{38,48} 6 
{38,63} 6 
{48,49} 6 
y {49,63} 6 



9 

128 

v5 

3 

3 

8V2 

16 

3 



Vi 
3 

9 

25 
256 

9 



1- 



2^ 



1 
V2 



26^/2+27tan-l(2^/2) 
277r 

^ (9 + 2V37r) 

52+27\/2sec-i(3) 
48 Ve 



9 

4 

3 27r 

1 

2 

V3 
2 

92+75sin-i(|) 
757r 



0.707107 \ 

0.785398 

0.825312 

0.828450 

0.842035 

0.888889 

0.769800 

0.608998 

0.500000 

0.866025 

0.685627 



20 

25 
48%/2 



^ (^5 + 3v^csc-i (^^)^ 0.886838 

92+75cos-i(|) 



16 
25 
24^2 

9V37r 
80 

3 
^/Io 

3 
5^2 



49 
384 



80^5 

14 
15 

7 
5\/5 

75+3 CSC 



37r 



1 



2^ 



^ ^ ( 7 

7^7 16 ' ' 



49 
72^/2 

27r 
7 



5 
6 

22%/6+49sin-i(2^) 

49^? 

- 4v^CSC"^ 

3(22'V6+49cos-i(f )) 

112^7 
V7+4csc-i(^) 

47r 



A)) »■ 



0.903076 
0.933333 
0.626099 

0.504975 
0.408248 
0.912871 

0.596820 
915544 
0.930129 
0.440596 
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We see that there are fewer rows in p6|) than in ()14|). obtained by the alternative form of 
partial transposition. 



1. Boundary states 

Our analysis of the US-lengths of the corresponding boundary states yielded: 

^{3,6} I I I 0.333333^ 

{4,18} i 0. 

8} 175 I ^ 0.942809 

{6,15} I i I 0.75 

{8,9} ^ ^ I 0.666667 

{8,22} ^ ^ f 0.666667 

{8,29} ^ ^ i 0.25 

{9,15} 75 I ^ 0.866025 

{15,22} ^ ^ I 0.2 

{15,29} ^ I 0.75 

{18,24} I i t 0.8 

{18,35} ^ i I 0.833333 

{24,25} ^ ^ I 0.8 

{24,46} ^ ^ t 0.8 

{24,57} ^ ^ I 0.166667 

{25,35} ^ I ^ 0.745356 

{35,46} ^ ^ f 0.142857 

{35,57} ^ 4^ I 0.833333 

{38,48} ^ 1 ^ 0.699854 

{38,63} f i I 0.875 

{48,49} ^ 4^ f 0.857143 

^{49,63} ^ 1 4 0.661438^ 
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C. Bi-PPT 



Here, we obtain the probabilities that an 8 x 8 density matrix will have a positive partial 
transpose, under both forms of partial transposition employed immediately above. Our 
results were: 



{1,20} 


28 


{3,6} 


12 


{3,20} 


8 


{3,36} 


4 


{6,8} 


2 


{6,15} 


2 


{8,9} 


2 


{8,20} 


4 


{8,22} 


2 


{8,31} 


2 


{8,36} 


4 


{8,42} 


2 


{8,49} 


4 


{8,53} 


2 


{9,15} 


2 


{15,22} 6 


{15,29} 8 


{18,24} 6 


{18,35} 4 



1 

4 

1 

3^2 
1 

4 

1 

3^2 
_9 

128 A/ 2 

v5 

3 
3 

3 

8 

3\/3 
16 

3^3 
16 

3 
8^2 

3V3 
16 

3 
8^2 
3 



6\/2 

VI 

3 
9 



TT 



4 
1 
V2 
5 
6 

5 

26^/2+27tan-l(2^/2) 
277r 

^ (9 + 2^7r) 

52+27\/2sec~^(3) 



48^6 



1 

2 

3 

7 
9 



25 
256 



TT 



9 

7 

3^ 
4 

3^3 

3 27r 

1 

2 

2 

92+75sin-i(|) 
75^ 



0.785398 
0.707107 
0.833333 
0.883883 
0.825312 

0.828450 
0.842035 
0.577350 
0.666667 
0.777778 
0.816497 
0.888889 
0.952579 
0.769800 
0.608998 
0.500000 
0.866025 
0.685627 



20 



^ (5 + 3V5CSC-1 (^^)) 0.886838 
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/ 



{20,35} 2 

{20,48} 4 

{22,48} 4 

{22,63} 10 

{24,25} 8 

{24,46} 2 

{24,57} 2 

{25,35} 4 

{31,35} 2 

{35,36} 4 

{35,38} 4 

{35,46} 2 

{35,57} 2 

{36,48} 4 

{42,48} 4 

{48,49} 2 

{48,51} 2 

{48,53} 4 

\{49,63} 8 



20 



49 



192 

192 
8 

25 
48^/2 

Vi 

16 
25 
24a/2 

9V37r 
80 

9v^7r 
80 

3 
5V2 

3 
5V2 

3 
^/Io 

3 
5^2 

49 
384 

49 
384 

49 

49 

72 

49 
72V2 

27r 
7 



TM (-25 + 24a/30) 
2^(l2 + 7V6sin-i 

(-9 + 7V42) 
^(7 + 4^csc-i 



4 

147 



92+75 cos-i(l) 
8OV5 



14 

15 

7 

a/5+3 CSC 



■(A) 



37r 



14t/5+27csc-i(^;^) 



277r 

5 



24 



1 

14^5+27 cos-i (I) 
24a/6 

I 



2^/6+7sin-l(2^) 
7^? 

22A/6+49sin-i(^2^) 
497r 

3(2V6+7cos-i(|)) 

iWt 



1 - 



3(22V6+49cos-i(|)) 
lliTT 

47r 



0.985680 ^ 
0.903278 

0.989529 
0.915544 

0.903076 

0.933333 
0.626099 

0.504975 
0.636783 
0.809819 
0.918793 
0.408248 
0.912871 

0.469522 

0.596820 
0.731739 
0.801610 
0.930129 
0.440596 J 
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1. Boundary states 



The lengths, separable lengths and separabihty probabihties of the corresponding (exte- 
rior/degenerate spectra) boundary states are given in the following array: 



/ 



{1,20} 

{3,6} 

{3,20} 

{3,36} 

{6,8} 

{6,15} 

{8,9} 

{8,20} 

{8,22} 

{8,31} 

{8,36} 

{8,42} 

{8,49} 

{8,53} 

{9,15} 

{15,22} 

{15,29} 

{18,24} 



1 

2 


U 


U 


3 


1 


]^ 


4 


4 


3 


1 


1 


1 


2 


4 


2 


3 


1 


1 


4 


4 


3 


3 


1 


2^2 


4^2 


2 


3 


1 


1 


3 


3 


4 


4 


3v^ 


V3 


2 


8 


4 


3 


V3 


V3 


1 


2 


8 


4 


3^3 
8 


n 





3^3 


V3 


1 


8 


8 


3 


3^3 


V3 


2 


8 


4 


3 




V3 


2 


8 


4 


3 


3V3 


V3 


2 


8 


4 


3 


V3 


V3 


1 


2 


8 


4 


1 


1 


V3 


x/3 


2 


2 


5 


1 


1 


2^6 


2^/6 


5 




A/f 


3 


V 3 


2 


4 


5 


1 


4 


8 


2 


5 


3 


1 


5 


10 


4 


6 



\ 



0.5 

0.333333 

0.942809 
0.75 

0.666667 



0.333333 

0.666667 
0.666667 
0.666667 
0.25 



0.75 
0.8 

0.833333 



(19) 



20 



{20,35} 
{20,48} 



24 
7 



{22,48} ^ 



3 
3 

3- /3 



{22,63} 
{24,25} 
{24,46} 
{24,57} 
{25,35} 
{31,35} 
{35,36} 
{35,38} 
{35,46} 
{35,57} 
{36,48} ^ 
{42,48} 4 
{48,49} 
{48,51} 
{48,53} 



^{49,63} ^ 



1 


6 


4 


7 






V 6 




2 


V 


1 


7 


4 


8 




4 


2 


5 




4 


2 


5 




1 


8 


6 


1 


vz 


2 


3 


1 


Vb 


2 


3 


vl 


1 


4 


6 




5 


4 


6 


V 5 


X 


4 


7 




5 


4 


6 


1 


2^6 


2 


7 


1 


2^6 


2 


7 










1 


8 


7 


\/2T 


6 


4 


7 


1 




2 


4 



^ 0.912871 
0.857143 
0.92582 
0.875 
0.8 
0.8 

0.166667 
0.745356 
0.745356 
0.166667 
0.833333 
0.142857 
0.833333 
0.699854 
0.699854 
0. 

0.142857 
0.857143 



D. Tri-ppt 



Now, we derive the probabilities that an 8 x 8 density matrix will have a positive partial 
transpose, not only under both forms of partial transposition previously employed, as in 
sec. IIVCI but also under a third (independent) form obtained, first, applying a certain 8x8 
permutation matrix (Q, eq. (3)]) to the original 8x8 density matrix, then transposing in 
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place the resultant four 4x4 blocks. We obtained the following results. 



{1,20} 


288 


{3,6} 


12 


{3,20} 


4 


{3,22} 


4 


{3,36} 


4 


{6,8} 


2 


{6, 15} 


2 


{8,9} 


2 


{8,20} 


4 


{8,22} 


2 


{8,31} 


2 


{8,36} 


6 


{8,42} 


2 


{8,51} 


2 


{8,53} 


2 


{9,15} 


2 


{15,22} 6 


{15,29} 10 


{18,24} 6 


{18,35} 4 



1 

4 




TT 
4 




0.785398 


1 

3\/2 




1 

V2 




0.707107 


1 
4 




5 
6 




0.833333 


1 

4 




2 
3 




0.666667 


1 

3\/2 

9 ; 
128 Y 




5 

4^2 

26\/2+27tan-l(2\/2) 
277r 




0.883883 
0.825312 


3 




^ (9 + 2v^7r) 




0.828450 


3 

8v^ 




52+27^/2 scc-l(3) 

48 Ve 




0.842035 


3 

8 




1 

V3 




0.577350 


3\/3 
16 




2 
3 




0.666667 


3^3 
16 




7 
9 




0.777778 


3 
8%/2 




pi 
V ^ 




0.816497 


3\/3 
16 




8 
9 




0.888889 


3 

8^2 




7 
3^6 




0.952579 


3 
8 




4 
3^3 




0.769800 


TT 

6a/2 
3 




1 1 V3 
3 27r 

1 

2 




0.608998 
0.500000 


2v^ 
9 




V3 
2 




0.866025 


25 / 




92+75 sm-l(l) 




0.685627 


256 Y 


757r 




20 




^ (^5 + 3V5csc-^ 


m 


0.886838 



22 



4Vl 



/ {20,35} 2 

{20,48} 6 -,92 

{22,63} 12 ^ 

{24,25} 10 

{24,46} 2 ^ 

{24,57} 2 ^ 

{25,35} 4 ^ 

{31,35} 2 

{31,48} 2 

{35,36} 2 ^ 

{35,38} 4 ^ 

{35,46} 2 

{35,51} 2 ^ 

{35,57} 2 

{36,48} 4 

{42,48} 4 

{48,49} 2 ^ 

{48,51} 2 ^ 

{48,53} 4 ^ 

\^{49,63} 8 



m (-25 + 24v^) 
^(l2 + 7V6sin- (^) 
3^(7 + 4v^csc-i 



92+75 cos- 



-l/'3^ 



80\/5 



9v^7r 
49 

~92 

3 



3_ 
3_ 

3 

5v^ 

49 
384 

49 
384 

49 

49 
2v^ 

49 



14 
15 

7 

5^5 
v%+3c 



37r 



14^5+27 csc-l(^^) 



277r 



4 

147 



(-9 + 7^42) 



24 



(v^ + 3cos-i(|)) 



1 

14^5+27 cos-l(|' 
2475 



2V6+7sin-i('2^1 



TT 



7lT 

sm 



22^6+49 2v:i'| 



TT 



2 
7 

2 " 497r 

3(2^/6+7cos-l(f )) 
16^7 

1 _ 

7 

3(22x/6+49cos-i(f )) 
112^7 



2n 
7 



V7+4C 



.7? J 



4n 



Here there are 
simply the "B 



only four generator pairs yielding the probability 
i-PPT" case fsec HvT^ . 



0.985680 \ 
0.903278 
0.915544 
0.903076 
0.933333 
0.626099 

0.504975 
0.636783 
0.989529 
0.809819 
0.728612 
0.408248 
0.918793 
0.912871 

0.469522 
0.596820 
0.731739 
0.801610 
0.930129 
0.440596 / 



|, while there were eight in 
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1. Boundary States 



Now, we obtained from the analysis of the one-dimensional (exterior) boundary states 
the results: 

/ {1,20} I 0. \ 

{3,6} III 0.333333 

{3,20} i 1 1 0.5 

{3,22} I 0. 

{3,36} III 0.333333 

{6, 8} ^ i V 0.942809 

{6,15} I I I 0.75 

{8,9} ^ ^ I 0.666667 

{8,20} f f i 0.25 

{8,22} ^ 0. 

V ^ (21) 

{8,31} ^ ^ I 0.333333 

{8,36} ^ ^ I 0.666667 

{8,42} ^ ^ I 0.666667 

{8,51} ^ ^ I 0.666667 

{8,53} f f i 0.25 

{9,15} ^ i ^ 0.866025 

{15,22} ^ ^ I 0.2 

{15,29} ^ I 0.75 

{18,24} I i I 0.8 

V{18,35}^ I I 0.833333/ 
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{20,35} 
{20,48} 
{22,63} 
{24,25} 
{24,46} 
{24,57} 
{25,35} 
{31,35} 
{31,48} 
{35,36} 
{35,38} 
{35,46} 
{35,51} 
{35, 57} 
{36,48} 
{42,48} 
{48,49} 
{48,51} 
{48,53} 
V{49,63} 



3 

5 

J_ 

24 

2 

kl 

8 



4 
3 
2^5 
3 

2v^ 

7_ 

12 

3%/f 



2 

Vi 

2 
7 
4x/6 
7 

4v^ 

Vi 

8 

Vi 

8 



2 

77 



3 
3 





1 


4 


6 










1 


4 


7 


VIE 


5 


4 


6 




5 


4 


6 


1 


2^ 


2 


7 


1 


2^6 


2 


7 








vT 


1 


8 


7 




6 


4 


7 


1 


77 


2 


4 



0.912871 

0.857143 
0.875 
0.8 
0.8 

0.166667 



^ 0.745356 
0.745356 



f 0.92582 
0.166667 
0. 

0.142857 
0.833333 
0.833333 



0. 

0.142857 
0.857143 
0.661438 
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2. Length of Triseparability-Nontriseparability Interior Boundary 



In the following array, we present the HS-length of the common (interior) border sepa- 
rating the triseparable states from the non-triseparable ones for each specific scenario. 



f U,20| 


1 


|3,6| 


3 
4 


|3,20| 


3 
4 


13 221 


1 

2 


|3,36| 


3 
4 


|6,8| 


3 
2 


|6, 15} 


5 
4 


{8,9} 


V3 
2 


{8,20} 


Q 
O 


{8,22} 


5^3 

Q 
O 


{8,31} 


5^3 
8 


{8,36} 


sfZ 

2 


{8,42} 


VI 

2 


{8,51} 


V3 
2 


{8,53} 


5^ 
8 


{9,15} 


1 


{15,22} 


3^ 
2 


{15,29} 


7 
2^6 


{18,24} 


3 
2 


i^{18,35} 


5 
4 



(22) 
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I {20,35} 


3 A / T 

V 10 


{20,48} 


5 
4 


{22,63} 


1 


{24,25} 


4 


{24,46} 


4 


{24,57} 


Q 
O 


{25,35} 


1 


{31,35} 


3 
2 


{31,48} 


2 


{35,36} 


4 


{35,38} 


3v^ 
2 


{35,46} 


11^ 

4 


{35,51} 




{35,57} 


4 


{36,48} 


1 


{42,48} 


3 

2 


{48,49} 


13^ 

8 


{48,51} 


13^ 
8 


{48,53} 


Ft 

V 3 


^{49,63} 


3 
4 



V. THE QUTRIT-QUTRIT CASE n = 9 

Here we only have — since 9 = 3^ — one option available for computing the partial 
transpose, that is transposing in place the nine 3x3 blocks of the 9x9 density matrices. 
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We obtained the results: 



/ {1,13} 

{3,11} 8 

{3,13} 8 

{3,25} 8 

{8,13} 8 

{8,25} 8 

{11,15} 4 

{11,24} 2 

{13,24} 2 

{13,35} 4 

{15,16} 4 

{15,44} 4 

{15,55} 4 

{16,24} 2 

{20,24} 2 

{24,25} 2 

^{24,27} 2 



360 



16 
81 

32^2 
243 

16 
81 

32\/2 
243 

_8_ 
27 

4V2 
27 

16\/27r 

243 

162 
2571 
162 

16 



15V15 

128^2 

729 
32^ 



81 
128 



243 
25 

324 V 2 

25 /5 
324 Y 2 

50^2 
243 

50V2 
243 



TT 



7r 
4 

1 

V2 
5 
6 

5 

4V2 
1 

%/3 

3 

1 I 3V3 

3 47r 



l + isin 



Ml) 
A(-2 + 5v/5) 



^ (^5 + 3a/5csc-i ^ 



VV5 



^(9y3 + 47r) 

n 

12 

n 

16 

4+5sin-i(|) 
5n 

92+75sm-i(|) 
757r 

1 - ^ 

^ 5^5 

92+75cos-i(|) 
80^5 



0.785398 \ 

0.707107 

0.833333 

0.883883 

0.577350 

0.816497 

0.746830 

0.863648 

0.979236 

0.886838 

0.879838 

0.916667 

0.687500 

0.549815 

0.685627 

0.821115 

0.903076 / 
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/ {24,44} 2 

{24,46} 2 
{24,55} 2 
{24,59} 2 
{24,70} 2 
{24,72} 2 

{25,35} 4 
{35,46} 4 
{35,70} 4 

{38,48} 8 
{38,63} 8 
{40,63} 8 
{48,49} 8 

{49,63} 8 

{53,63} 8 

{63,64} 8 

\{63,66} 8 



100V2 


1 


243 


V5 




11 


81 


15 


50V2 


2 


243 


V5 


lOx/lO 


14 


81 


15 


50-\/2 


11 


z4to 


ov 5 


100\/2 


7 


243 


5a/5 


47r 


V5+3 


15^3 




32\/f 


1 


27 


VE 


32^2 




135 


V 6 



49 
486 

512 
567^ 
512 

567v^ 

196^2 
729 

1287r 
567 

1287r 
567 

512V2 
1701 

512\/2 
1701 



37r 



TT 



22V6+49sin-l(^^j 
49^;^ 



(7 + 4v/7c=c-'(^)) 



49 I \/l4 
192 ~'~ 3 
3(22^6+49 cos-l(l)) 

U277 
V7+4csc-i(-^) 

47r 

19-^/7+48 csc-i(^) 

487r 

7 



64 



19^7+48 cos-i (I) 
6475 



0.447214 \ 
0.733333 
0.894427 
0.933333 
0.983870 
0.626099 

0.504975 
0.408248 
0.912871 

0.596820 
0.915544 
0.992011 
0.930129 

0.440596 
0.563412 
0.795378 
0.938690 / 
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A. Boundary States 



Here, we have for the one-dimensional sets of (exterior) boundary states the results 



/ {1,13} 


4 
9 








0. 


13, 111 


2 
3 


2 
9 


1 
3 


0.333333 


{3, 13} 


4 
9 


2 
9 


1 
2 


0.5 


{3,25} 


2 
3 


2 
9 


1 
3 


0.333333 


{8,13} 


4 
3\/3 


1 

3v3 


1 

4 


0.25 


{8,25} 


1 
Vo 


2 

o V "J 


2 
3 


0.666667 


{11,15} 


8 

9i/3 


4 

9 


^/3 
2 


0.866025 


{11,24} 


5 
18 


2 
9 


4 

5 


0.8 


{13,24} 


5 
9 


2^5 
9 


2 

V5 


0.894427 


{13,35} 


4 
15 


2 
9 


5 
6 


0.833333 


{15,16} 


9 


3 


3 
4 


0.75 


{15,44} 


9 




3 
4 


0.75 


{15,55} 


10^ 
9 


9 


1 
5 


0.2 


{16,24} 


5 
9 


4 
9 


4 

5 


0.8 


{20, 24} 


5 
9 


4 
9 


4 

5 


0.8 


{24,25} 


9 


9 


1 

5 


0.2 


V{24,27} 


9 


2x/l0 
9 


4 
5 


0.8 
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/ {24,44} 4^ ^ I 0.166667^ 

{24,46} ^ ^ I 0.2 

{24,55} ^ ^ I 0.8 

{24,59} ^ ^ I 0.8 

{24,70} ^ ^ I 0.8 

{24,72} ^ ^ I 0.166667 

{25,35} ^ I f 0.745356 

{35,46} ^ ^ 1 0.142857 

{35,70} ^ ^ I 0.833333 ■ 

{38,48} ^1 ^ 0.699854 

{38,63} i I I 0.875 

{40,63} § ^ ^ 0.935414 

{48,49} ^ ^ f 0.857143 

{49,63} ^ I ^ 0.661438 

{53,63} ^ I ^ 0.661438 

{63,64} ^ ^ I 0.125 

V{63,66} ^ ^ I 0.875 / 

B. Length of Separability-Nonseparability Interior Boundary 

In the following arrays, we present the HS-length (and now a numerical approximation 
to it) of the common border separating the states lacking a PPR from those that possess a 
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PPT for each specific scenario. 



' {1,13} 


4 

9 


0.444444 ^ 


r r» 111 

{3,11} 


4 
973 


0.2566 


{3,13} 


A 

9 


0.444444 


{3,25} 


4 

973 


0.2566 


{8,13} 


^ (-2V3 + 3V6) 


0.287731 


{8,25} 


(1173 + 12^33) 


0.296254 


{11,15} 


4 

9 


0.444444 


{11,24} 


4 
11 


0.363636 


{13,24} 


4 

Ti 


0.363636 


{13,35} 


y 

18 


0.388889 


{15,16} 


(23\/6 + 12 V138) 


0.317724 


{15,44} 


^ (3 + v^) 


0.403666 


{15,55} 


^(-5^6 + 8^) 


0.290017 


{16,24} 


4 
9 


0.444444 


{20,24} 


4 
9 


0.444444 


{24,25} 


117VlO+40\/390 


0.330462 


3510 


U24,27} 


117VT0+40V390 
3510 


0.330462 y 
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( {24,44} 
{24,46} 
{24,55} 
{24,59} 
{24, 70} 
{24, 72} 
{25,35} 
{35,46} 
{35,70} 
{38,48} 
{38,63} 
{40, 63} 
{48,49} 
{49,63} 
{53,63} 
{63,64} 
\^{63,66} 



4 
■ 135 



(3\/T0-5 
^(4 + v^) 



15) 



18 

117\/10+40a/390 

3510 

i(4 + v^) 

117^/lO+40V390 
3510 



4 
135 



(3v^ - 5v^) 



-Ti9(7v^-12V2T) 

4(59yT5+15\/885) 
7965 

4 
9 

12 
29 

12 

29 

415y2T+84\/l743 



15687 

4 
9 

4 
9 

333V7+56V777 
6993 

333^/7+56V777 
6993 



0.292684 
0.397904 
0.330462 
0.397904 
0.330462 
0.292684 
0.444444 
0.295027 
0.338853 
0.444444 
0.413793 
0.413793 
0.344789 
0.444444 
0.444444 
0.349209 
0.349209 



VI. THE CASE n = 10 



A. 5x2 Decomposition 



We first compute the partial transpose by transposing in place the (2^) four 5x5 blocks 
of our set of two-dimensional 10 x 10 density matrices. Our analysis yielded for the HS-total 
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volumes and PPT-probabilities 



/ ri on"! 

|1,29| 


480 


4 

25 


TT 

4 


{3,27} 


4 


T5~ 


I 

73 


{3,29} 


12 


4 
25 


5 
6 


{3,49} 


1 O 


8\/2 
75 


5 
4a/2 


|«, zy| 


A 


6 

25 


1 

V3 


{8,31} 


o 

8 


25 


7 

9 


18 491 


4 


3^2 

25 


V 3 


|», 00 


A 


3V3 
25 


8 
9 


{8,64} 


Q 
O 


3\/2 
25 


7 

3\/6 


18 68\ 


4 


6 
25 


4 
3^3 


|iO, dij 





32^ 
75 


1 
2 


i 1 ^ 331- 




^Vl 

25 


3 
4 


-fl5 64]- 


6 


32^2 
225 


V3 
2 


{15,72} 


o 
z 


25 


11 
12 


{15,81} 


6 


32^2 
225 


9a/3 
16 


{15,87} 


2 


32^ 
75 


11 
16 


{24,33} 


8 


V2 
3 


1 

V5 


{24,81} 


8 


1 

3^2 


2 

x/5 


{27,35} 


8 


125 


14^+27 CSC" 
277r 


^{27,48} 


2 


49v^ 
300 


^ (|l2 + 7\/6sin 



0.785398 \ 
0.707107 
0.833333 
0.883883 
0.577350 
0.777778 
0.816497 
0.888889 
0.952579 
0.769800 
0.500000 
0.750000 
0.866025 
0.916667 
0.974279 
0.687500 
0.447214 
0.894427 

0.636783 
0.903278 y 
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^ {29,48} 


6 


49VT 
300 


^ (-9 + 7^/42) 


0.989529 


\ 


{29,63} 


4 


128 
175>/7 


^(7 + 4v/7csc-i U)^ 


) 0.915544 




{31,63} 


8 


128 
175a/7 


49 _|_ \/l4 
192 3 

4 , 

9 V2 


0.992011 




131 80\ 


6 


243 
800 


925046 




{33, 80} 


6 


243 

oUU 


^f, (-8 + 27a/2) 


0.993704 




{33,99} 


8 


27 


|(3 + 5sm-i (!)) 


0.932626 




{35,36} 


8 


24^/2 
125 


/ — 1 / \ 

14x/5+27cos-i(|) 

24^6 


0.918793 




{36,48} 


2 


49 /t"- 
600 Y 2'' 


2V6+7sm-i(2^) 
7n 


0.469522 




{40,48} 




600 Y 2 


22^6+49 sin-i(^^) 






6 


497r 


0.596820 




{48,49} 


6 


49 
225 


1 '^V'^ 
^ 7 


0.801610 




{48,51} 


6 


49V2 
225 


3(22^6+49 cos-l(l)) 
112\/7 


0.930129 




{49,63} 


4 


327r 
175 


V7+4csc-i(^) 
47r 


0.440596 




{55,63} 


4 


327r 
175 


1 Q. /7_i_/ie ^c^ — 1 1 4 \ 
lyy /-r4cScsc 1 —7= 1 

\ V ' / 

487r 


0.563412 




{63, 64} 


8 


128^/2 
525 


1 64 


0.795378 




{63,68} 


4 


128a/2 
525 


19^7+48 008-1(1) 
64 a/2 


0.938690 




{64,80} 




2437r 


4^/2+9sin-l(i^) 






6 




0.416417 




8OOV2 


97r 




172 801 


2 


2437r 


ll^+sin-(i#) 


534977 




800 V2 


TT 




{80,81} 


6 


27 
50V2 


I 4^2 


0.790487 




{80, 87} 


2 


27 
50V2 




0.945383 




l^{81,99} 


8 


V5tt 
9 


3+5 sin-i(l) 
57r 


0.395819 


/ 
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B. 2x5 Decomposition 

Now, we compute the partial transpose by transposing in place the (5^) twenty- five 4 
blocks of the 10 x 10 density matrices. 



^{3,6} 


16 


8\/2 
75 


1 
V2 


0.707107 


{4, 18} 


480 


4 
25 


4 


0.785398 


{6,8} 


2 


200 Y 2 


26^2+27 tan-l(2V2) 
277r 


0.825312 


{6,15} 


2 


75 


^ (9 + 2V37r) 


0.828450 


{8,9} 


2 


3^2 

25 


52+27A/2sec-i(3) 
48 a/6 


0.842035 


{8,22} 


6 


3^3 
25 


8 
9 


0.888889 


{8,29} 


6 


6 
25 


4 
3a/3 


0.769800 


{9,15} 


2 


4\/27r 
75 


1 1 V3 
3 ^ 27r 


0.608998 


{15,22} 


6 


75 


1 
2 


0.500000 


{15,29} 


6 


32 V2 

225 


V3 
2 


0.866025 


{18,24} 


4 


-a/^vt 

16 Y 2 


92+75 sin-i(l) 
757r 


0.685627 


{18,35} 


4 


36vf 
125 


^(5 + 3V5csc-i (^)^ 


) 0.886838 


{24,25} 


4 


1 


92+75 cos-i(l) 


0.903076 


3^2 


80^/5 


{24,46} 


4 


1 

Vw 


14 
15 


0.933333 


V{24,57} 


4 


V2 
3 


7 

5a/5 


0.626099 



36 



{25,35} 4 

{35,46} 4 

{35,57} 4 

{38,48} 6 

{38,63} 6 

{48,49} 6 

{48,78} 2 

{48,93} 2 

{49,63} 6 

{63,78} 2 

{63,93} 2 

{66,80} 8 

{66,99} 8 

{80,81} 8 

\^{81,99} 8 



9V3n 
125 

25 

24^2 
125 



600 Y 2 
128 



V5+3csc-i(^) 



175^7 

49^2 
225 

!Vl 

25 

98 
75V3 

32tt 
175 

128v^ 
75 

128\/2 
525 

2437r 
800 V2 
4V5 

27 

27 



50V2 

■v/Stt 
9 



1 

5 
6 

22^/6+49sin-l(^2ySj 



497r 



^7 

16 ' ' 



4\/7csc-^ 

3(22^6+49 cos-i(f)) 



112a/7 



20 
21 

10 

%/7+4csc" 



47r 



1 

2a/2 



2_ 

2 

172\/2 
243 



4 

V7 



A(3 + 5sin-i(|)) 

3+5sin-i(|) 
57r 



0.504975 
0.408248 
0.912871 

0.596820 
0.915544 
0.930129 
0.952381 
0.539949 

0.440596 
0.353553 
0.935414 

0.534977 
0.932626 
0.945383 
0.395819 
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C. Bi-PPT 

Now, we require that the 10 x 10 density matrices be positive under both the forms of 
partial transposition employed immediately above. We obtained the (extensive) results 



( {1,29} 
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7r 
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0.785398 \ 


{3,6} 
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1 

V2 




0.707107 


{3,29} 
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25 


5 
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0.833333 


{3,51} 


6 


8\/2 
75 


5 
4V2 




0.883883 


{6,8} 


2 


200 V 2 


26^2+27 tan- 
277r 


-i(2V2) 


0.825312 


{6,15} 


2 


75 


^ (9 + 2V37r) 


0.828450 


{8,9} 


2 


3\/2 
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52+27 V2sec- 
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0.842035 


{8,22} 
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0.888889 
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1 
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0.577350 
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0.777778 
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0.666667 
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3\/3 




0.769800 
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0.952579 
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1 I 

3 27r 




0.608998 


{15,22} 10 
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2 




0.500000 


{15,29} 12 
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V3 
2 




0.866025 


{15,33} 6 


25 


3 
4 




0.750000 


{15,72} 2 


25 


11 

12 




0.916667 


V {15,81} 4 


32^2 
225 


9\/3 
16 




0.974279 ) 



(28) 
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{15,87} 2 

{18,24} 4 

{18,35} 4 

{24,25} 4 

{24,33} 8 

{24,46} 4 

{24,81} 8 

{24,89} 2 

{25,35} 4 

{27,35} 6 

{27,48} 2 

{29,48} 6 

{29,63} 10 

{31,63} 6 

{31,80} 6 

{33,80} 6 

{33,99} 16 

{35,36} 8 

{35,46} 4 

V{35,57} 4 
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1 /5 



125 
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3^2 
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VlO 
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3\/2 
V2 
3 

9\/37r 
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9\/3; 
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300 
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175v^ 

243 
800 
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4-^/5 
27 

24^2 
125 

!Vi 

25 

24x/2 
125 
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16 

92+75sin-i(|) 

757r 

ll(5 + 3V5csc-i (^)) 

92+75cos-i(|) 
80V5 

1 

14 
15 

2 

v/5 
7 
5^5 

V5+3csc-i(^) 



37r 



14^5+27 csc-i(^;^) 
277? 

^(l2 + 7V6sin- (^)) 

Ti7 (-9 + 7V42) 

1^(7 + 4a/7csc- (^)) 



49 _|_ a/T4 
"l92 ~'~ 3 
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V2 

213 (-8 + 27V2) 
(3 + 5sin-^(|)) 
^(1) 



_3^ 
20 

14a/5+27 cos 
24^6 

1 

V6 



0.687500 

0.685627 

0.886838 

0.903076 
0.447214 

0.933333 
0.894427 
0.626099 

0.504975 

0.636783 
0.903278 

0.989529 
0.915544 

0.992011 

0.925046 
0.993704 
0.932626 
0.918793 
0.408248 
0.912871 
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{36 



{38 
{48 
{48 
{48 
{48 
{48 
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{55 
{63 
{63 
{63 
{63 
{64 
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{80 
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49} 2 
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225 

49V2 
225 



25 

98 
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525 
128 
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800v^ 
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21 
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2V2 

2 

4V2+9sin-i(4^) 
9^? 

il#+sin-(l^] 



0.469522 

0.596820 
0.731739 
0.930129 
0.801610 
0.952381 
0.539949 

0.440596 

0.563412 

0.795378 ■ 

0.938690 

0.353553 

0.935414 

0.416417 
0.534977 



^ (4^2 + 9 cos-i (I)) 0.735870 
1 - ^ 0.790487 

5472^16^"'' {9) 



3+5 sin 



0.945383 
0.395819 



The probability 0.993704 (corrresponding to the pair of Lie generators numbered {33, 80}) 
is the largest of any recorded in all our results above. (This also occurs in (j2SI)-) 



VII. ANALYSES OF SCENARIOS WITH MORE THAN TWO PARAMETERS 

We found it considerably simpler to extend the Jakobczyk-Siennicki model ^ from two- 
qubit systems {n = 4) to higher n — as illustrated above — than to extend it from two- 
dimensional sections of Bloch vectors to m-dimensional sections (m > 3), even just for the 
case n = 4. (However, we were able to obtain a highly interesting set of exact HS separability 
probabilites for certain m = 3, n = 4 systems, using the Jaynes maximum-entropy principle 
in conjunction with the new integration over implicitly defined regions feature of Mathe- 
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matica, in 0, Fig. 11].) A fortiori, it appears that the determination of the HS separable 
volume of the /i/i^een-dimensional convex set of 4 x 4 density matrices — conjectured on the 
basis of an extensive quasi- Monte Carlo a„aly=i= in to be (iVVS^ « 27 3707 • 10- 
eq. (41)] — would have to proceed in some quite different analytic fashion to that pursued 
here. (Based on our experience in the above-reported analyses, it appears to be a necessary 
condition for obtaining exact HS separability/PPT-probabilities that explicit formulas be 
available for the eigenvectors of both the class of density matrices under consideration and 
of their partial transposes.) 



A. m = 3, n = 4 



We have been able to find, up to this point in time, that for the i/iree-dimensional two- 
qubit {m = 3,n = 4) scenarios generated by the four triads of Gell-Mann matrices {1, 4, 6}, 
{1,5,7} {2,4,7} and {2,5,6}, the volume of separable states is — having to resort to nu- 
merical methods — 0.478512 and of all the (separable and nonseparable/entangled) states, 
0.61685, for a separability probability of 0.775734. For the scenario {10,12,13}, the sep- 
arable volume remains the same, but the total volume is exactly | ^ 0.523599 for an HS 
separability probability of 0.913891. 



B. Two-Dimensional Boundaries of m = 3, n = 4 Systems 

Of course, it we restrict attention to the generic boundary states of the three-dimensional 
scenarios, we only have to perform two-dimensional computations. Thus, we were able to 
find that for the triadic scenarios {1, 3, 6}, {1, 3, 7}, {1, 3, 9}, {1, 3, 10}, {1, 4, 9} and {1, 5, 10}, 
amongst others, the HS-area of the states with degenerate spectra is | and that of the 
separable component of this area, one-half that value. For the triadic scenarios {1, 4, 6} and 
{1,5,7}, the separable component of the boundary states has area ^ and the total area is 
I (^V5+6 sin"^(-^) j for a separability probability of 0.165025. Also, for several scenarios (for 
instance, {3,4,9}), we have a total area of ^y^, a separable area of |, giving a separability 
probability of ~ 0.353553. 

Now, we present all our results of this type (two-dimensional exterior boundaries of three- 
dimensional scenarios) in the following array. The first column gives the corresponding 
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triad of Gell-Mann matrices, the second column shows the total HS-area of the boundary 
states, the third column gives the exact separability probability, and the last, a numerical 
approximation to the probability. (There may exist additional nontrivial scenarios, as we 
were not readily able to fully analyze all 2730 = 13 ■ 14 • 15 possible triads of 4 x 4 Gell-Mann 
matrices.) 

/ {1,3,6} f I 0.500000 \ 



0.165025 
0.616044 



{1'4'6} i(v/5 + 6sin-(^)) 

{3,4,6} |(v/5 + 6sm-(^)) H .^^^ 

{3,4,9} ^ ^ 0.353553 

{3,6,7} i(v^ + 6sin-i(^)) 1 



0.500000 



, VT5-87r+8tan-i , ., 

{3,6,8} ^ ^ 0.636114 



(29) 



{6,9,15} f .-2V,^s.^ucs.-^iVE) Q 207232 

{8,9,10} ^ (4 + v/7 + 2. + 8 cot'^ {V7)) ,^^JXTl-^{Vi) O-^^OOl^ 
V {9, 11,13} ^ ^ 0.0833333 / 

C. 9- and 15-Parameter Analyses (m = 9, 15, n = 4) 

Now, we sought to make some progress in obtaining the (conjecturally exact) Hilbert- 
Schmidt volume of the separable 4x4 density matrices, in both the 9-dimensional case of real 
density matrices and the 15-dimensional case of (fully general) cornplex density matrices. In 
both cases, we dispensed with the Bloch vector parameterization used in the above 

analyses (neither did we employ the integration over implicitly defined regions capabilities 
of Mathematica version 5.1), and adopted a simple, naive parameterization, in which the 
four diagonal elements of the density matrices were denoted a, 6, c, 1 — a — 6 — c and the off- 
diagonal (upper triangular) elements, aij+i(3ij (in the real m = 9 case, of course, all /3's equal 
zero). (In order to compare our results here with the HS- volume formulas of Zyczkowski and 
Sommers Q], the volumes we do report below are our initial volumes multiplied by factors 
of 2^ in the complex case, and 2^ in the real case. In all our earlier analyses above, we have 
simply taken the HS- volume element to equal 1.) 

In both (m = 9, 15) of these cases, we pursued the same analytical strategy. We required 
that the six principal 2x2 minors of the density matrices and/or their partial transposes 
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have nonnegative determinants. This is (only) one of the requirements for nonnegative- 
definiteness (cf. [l^ eq. (12)]). Ideally, we would also have required that the leading 
principal 3x3 minor have nonnegative determinant and also that the determinant of the 
matrix be nonnegative. But these last two requirements were too computationally onerous 
to impose (at least in our first round of efforts). So, our analytical strategy should yield 
upper bounds on the Hilbert-Schmidt volumes in these cases. 

1. 9- Dimensional Real Case 

When we only required that the six principal 2x2 minors have nonnegative determinants, 
we obtained for the volume the result ~ 0.00881215. (We can reduce [improve] this 
to ^^gg^ ^ 0.00712396 by modifying [ narrowing], to begin with, the integration limits 
over a single off-diagonal variable, so that in addition, to its corresponding 2x2 minor, a 
corresponding 3x3 minor also has a nonnegative determinant. If we narrow similarly a 

4 

second set of integration limits, this is further reduced to ~ 0.00362385. An attempt 
to add a third set of similar integration limits — corresponding to a 3 x 3 minor — did 
not succeed computationally.) Applying formula (7.7) of the Zyczkowski-Sommers study 



j5|, we obtain for the HS- volume of the 4x4 real density matrices, ~ 0.0016106. 

This is |j ~ 0.18277 times sma//er than our first, principal calculation {^^), so we have a 
considerable overestimation. 

If we additionally imposed the condition that the six principal minors of the partial 
transpose also have nonnegative determinants (only two of them being actually different 
from the original six), the result was ^ 0.00548249. (Note that 1120 = 2^ • 5 • 7 and 
99225 = 3'^ ■ 5^ ■ 7^.) So (taking the ratio) of this to ^f^, we obtain a crude estimate of the 
HS-separability probability of the real density matrices is 0.622151. 

Unfortunately, our upper bound (0.00548249) on the HS-volume of the separable real 
two-qubit states is larger than the (known) HS-volume (0.0016106) of the (separable and 
nonseparable) two-qubit states, so we have not yet succeeded in deriving a nontrivial upper 
bound on the separable volume. (The same will be the case in the immediate next analysis.) 




.4 
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2. 15-Dim.ensional Complex Case 



Now, when we again required that the six principal 2x2 minors have nonnegative deter- 
minants, we obtained for the corresponding Hilbert-Schmidt volume ygg^^g^g ~ 0.000121959. 
(Note that 7882875 = 32-53-72-11-13.) Formula (4.5) of fl gives us for the HS-volume of the 
15-dimensional convex set of two-qubit density matrices, the value §5^30500 ~ 1.12925 ■ 10~ . 
(The ratio of these two volumes — the measure of our overestimation — is ~ 10.7994.) 
Imposing (just as we did in the real 9-dimensional case) the further requirements that the 
six 2x2 minors of the partial transpose all have nonnegative determinants, we obtain a 
HS-volume of 304^375 ~ 0.0000620378. (Observe that 30435780375 = 3^ ■ 5^ ■ 7^ ■ II2 ■ 13^.) 
So, our crude separability probability estimate (less than in the 9-dimensional real case as 
conforms with our intuition) is ~ 0.508677. Based on certain numerical and theoretical 
considerations, the actual value of this (15-dimensional) separability probability has been 
conjectured to be 0, eq. (43)] 

22. 3 -72 -11 -1373 
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VIII. CONCLUDING REMARKS 



0.242379. (30) 



We have found the newly-introduced capability of Mathematica (version 5.1) for integra- 
tion over implicitly defined regions particularly useful for obtaining a very wide variety of 
separability (and positive-PPT) probabilities, particularly for low-dimensional (m = 2, 3) 
cases, essentially independently of the sizes (n) of the corresponding n x n density matrices 
analyzed. The use of such methods for cases m >= 4 appears, however — such as the 
two-qubit {n = 4) scenarios for the real (m = 9) and complex (m = 15) cases — to be 
particularly challenging. 

Eggeling and Werner studied the separability properties in a five-dimensional set of 
states of quantum systems composed of three subsystems of equal but arbitrary finite Hilbert 
space dimension. They are the states that commute with unitaries of the form U ®U ®U . 
In [1^, we evaluated the probabilities of an Eggeling- Werner state being biseparable, trisep- 
arable or having a positive partial transpose with respect to certain partitions. However, 
the Hilbert-Schmidt measure was not employed, but rather the Bures one jl5| . 
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